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Abstract. We analyze a time sequence of Inter-Network (IN) magnetograms observed at the solar disk center. Speckle recon- 
struction techniques provide a good spatial resolution (Of! 5 cutoff frequency) yet maintaining a fair sensitivity (some 20 G). 
Patches with signal above noise cover 60 % of the observed area, most of which corresponds to intergranular lanes. The large 
surface covered by signal renders a mean unsigned magnetic flux density between 17 G and 21 G (1 G= 1 Mx cm -2 ). The dif- 
ference depends on the spectral line used to generate the magnetograms (Fei /16302.5 A or Fei /(6301.5 A). Such systematic 
difference can be understood if the magnetic structures producing the polarization have intrinsic field strengths exceeding 1 kG, 
and consequently, occupying only a very small fraction of the surface (some 2%). We observe both, magnetic signals changing 
in time scales smaller than 1 min, and a persistent pattern lasting longer than the duration of the sequence (17 min). The pattern 
resembles a network with a spatial scale between 5 and 10 arcsec, which we identify as the mesogranulation. The strong depen- 
dence of the polarization signals on spatial resolution and sensitivity suggests that much quiet Sun magnetic flux still remains 
undetected. 
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1. Introduction 

Most of the solar surface appears as non-magnetic in tradi- 
tional magnetic field determinati ons (e. g., in the full-disk Kitt 
Peak magnetograms: LFones et alJll992l) . This so-called quiet 
Sun does not produce enough polarization to show up in the 
measurements. However, such lack of detection does not im- 
ply the non-existence or irrelevance of the quiet Sun mag- 
netism. On the contrary, the limited sensitivity of the stan- 
dard magnetograms and the large area covered by the quiet 
Sun point out that traditional measurements may easily over- 
look a large fraction of the solar magnetic flux. With various 
flavors and shades, this argument has been put forward many 
times during the past fifty years (e.g.lU nno 1959; Stenflo 1982; 
IZiridll987tlYi et alJll993t ISanchez Almeidalll998l 120031) . If 
the conjecture were correct and the quiet Sun carries a size- 
able fraction of the existing magnetic flux, then weak polar- 
ization signals should appear upon improvement of the sen- 
sitivity of the magnetograms. Such weak signals are actu- 
ally observed thanks to the last generation of solar spectro- 
polarimeters. When the noise is in the few G level and the angu- 
lar resolution abou t 1", then most of the solar surface becomes 
magnetic (e.g..|Grossmann-Doerth et alll99fJ:lLin & Rimmele 
ll999llLitesll2002l) . In addition, Hanle depolarization measure- 
ments of chromospheric lines also indicate the presence of an 
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ubiquitous magnetic field (e.g., tF aurober t-Scholl et alJ Il995t 
iBianda et alll999tlShchukina & Truiillo BuenofcOOsIT 

Little is known about the physical properties of these 
fields (magnetic flux, distribution of field strengths, structure, 
degree of concentration and tangling, relationship with the 
fields in the chromosphere and corona, etc.). The observational 
studies are still exploratory, aiming at setting up the scene. 
With this general purpose, we have observed the magnetic 
fields of an Inter-Network (IN) region with good angular 
resolution, yet maintaining a fair sensitivity. The study refers 
to low flux features in the interior of a network cell, excluding 
network patc hes for which an extensive literature exists 
(e.g., ISolankilfl993t IStenflolll994l IN magnetic fields were 
already discovered in the seventies dLivingston & Harvevl 
1 1 97.4 ISmithsonlll975l) . During the last decade, with the new 
observational and diagnostic capabilities, these IN fields have 
received increasing attention. So far only a fraction of the 
IN field has been detected. This conclusion follows from 
the presence of unresolved m ixed polarities in 1" angular 
resolution observations (e.g.. [Sanc hez Al meida et ail [l996; 
ISanchez Almeida & Litesl T2000I: lLitesll2002l) . The mixing of 
polarities reduces the polarization, making the IN magnetic 
structures difficult to identify. It is therefore to be expected 
that more locations with magnetic fields will be detected 
when increasing the spatial resolution. Some of the recent 
observational studies find IN fields having magnetic field 
strengths substantially lower than 1 kG ( Kelle r et al.l 1994; 
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1996; 



2000 



Linl 1 19951 iLin & Rimrneld Il999t iBianda et all 1 19981 Il999t 
Khomenko et alJ 120031) . On the contrary, different observ- 
ing and interpretation techniques indi cate the existence 
of stro ng kG magne ti c field s (jGros^manri^Doertiie^al. 

Sigwarth et al. 1999; Sanchez Almeida & Lites 
ISocas-Navarro & Sanchez Almeid a 2002). The in- 
consistency could be cured if the IN regions present 
a continuous distribution of field strengths. Depending 
on the specificities of the diagnostic technique, one is 
mostly se nsitive to a p a rticular part of such d i stribu- 
tion (see ICattaneol Il999t ISanchez Almeida & Litesl |2000; 
ISocas-Navarro & Sanchez Almeida! 12003). Lifetimes of IN 
feat ures have been measured to be between 0.2 to 7.5 hours 
JZhang et alJ Il998l) . Improving of the angular resolution 
and cadence redu ces the lifetimes, w hich come down to 
a few minutes jLin & Rimmeld Il999h . Such combination 
of short lifetimes with the large flux content makes IN 
regions a very efficient system to process magnetic fields, 
orders of magnitude more effective than both active re- 
gions (upliftin g 7 x 10 21 Mx day -1 during the maximum 
of the cycle, Harvev-Ang ld Il993l) an d ephemeral regions 
(~ 5 x 10 23 Mx day-'. lHagenaaibOOlh . 

Here we measure several basic physical properties of the 
IN fields. The novel side of the study lies in the good spa- 
tial resolution of the time series of magnetograms on which 
the measurements are based (075). The observations are pre- 
sented in Sect. |3 Then Sect. [3] describes how the high spatial 
resolution is obtained by speckle reconstruction, and how the 
magnetograms are produced from the reconstructed spectro- 
polarimetric maps. We analyze in detail the best snapshot of 
the sequence to characterize the properties of the IN fields 
(Sect.0}. The complete series allows us to study the time evolu- 
tion of the magnetic signals (Sect.[5j. In order to assess the reli- 
ability of our conclusions, we analyze in detail the consistency 
of the results, both internally and when compared with previ- 
ous measurements (Sect.[6}. Some o f the results p resented here 
were already advanced in a letter bv lDomfnguez Cerdefia et al] 
(200i. 



2. Observations 

The observations were obtained with the Gottingen Fabry- 
Perot Interferometer (FPI), which is a post-focus instrument of 
the German Vacuum Tower Telescope of the Observatorio del 
Teide (Tenerife, Spain). The target was a very quiet Sun region 
near disk center, which was selected using G-band video im- 
ages to avoid contamination from network magnetic concen- 
trations 1 . The optica l system derives fro m the early work by 
iBendlin et all i ll 9921) . iBendlinl (1 19931) . and lBendhn & Volkmet| 
( 1995). Its present setup is described in iKoschinskv et al] 
(2001). It includes two CCD cameras, operating simultane- 
ously, with an exposure time of 30 ms. Both CCDs have 
384x286 pixels, each pixel being 071 x 071 on the Sun, which 
corresponds to half the diffraction limit of the telescope at the 
working wavelength. The broad-band camera (CCD1) images 



1 Further evidence showing that the data belong to an IN region will 
be presented below. 



the field-of-view (FOV) through an interference filter (band- 
pass ~ 100 A, centered at ~ 6300 A). It provides the so-called 
speckle image, needed for the speckle reconstruction (Sect. 
13. U . The second camera (CCD2) gathers narrow-band images 
whose wavelength is selected by the FPI. This second beam in- 
cludes a Stokes V analyzer, which separates the left and right 
circularly polarized components of the light into different ar- 
eas of the CCD2 (for d etails, see Fig. 1 in IKoschinskv et all 
200lUKoschinskvl200lh . The FPI allows to scan in wavelength, 
taking several images per wavelength position. We select a 
spectral region around /I6302 A, which contains the iron lines 
Fe i ,16301.5 A (Lande factor g L = 1.67), Fei ,16302.5 A (g L = 
2.5), as well as the telluric line 2 ,16302.8 A. The FPI wave- 
length step was set to 31.8 mA, which suffices to sample the 
telluric line (5 wavelengths, with 2 images per position), and 
the two Fei lines (14 positions each, with 5 images per posi- 
tion). The dots in Fig^hidicate the exact wavelengths. A full 
scan with the properties described above renders 150 images, 
which is the maximum number allowed by the acquisition sys- 
tem. Each scan lasts 35 s, plus 15 s needed for storage onto 
hard disk. 

The FPI was adjusted to yield a band-pass of 
44 mA FWHM (Full Width Half Maximum). This finite 
wavelength resolution modifies the spectrum. Figure ^ shows 
the mean Stokes / profile in one of our flatfields, together with 
the spectrum of the disk center as publ ished in a stan dard solar 
atlas (Brault & Neckel 1987, quoted in lNeckelll 999). Note the 
good agreement between our spectrum and the atlas smeared 
with the theoretical band-pass of the FPI. 

The data used in the work belong to a short time series of 20 
scans taken in April 29, 2002. The cadence was some 50 s, so 
that the series spans 17 min. The seeing conditions were partic- 
ularly good, with a Fried parameter deduced from the speckle 
reconstruction between 13 and 14 cm. 



3. Data analysis 

3.1. Speckle reconstruction of the narrow-band 
images 

The reconstruction of the FPI narrow-band images is a two- 
step process, which begins by reconstructing the broad-band 
images. Using the 150 broad-band images taken in CCD1, 
one is able to obtain a reconstructed image with a spatial res- 
olution close to ff.'25 . The Gotting en code for speckl e re- 
construction is used as described in lde Boer et al] ill 9921) and 
deBoerl lll996). We use the spectral ratio method (Ivon der Luhel 
19841) to derive the seeing conditions and the amplitude cor- 
rection factors. The speck le masking method ( Wei geltl| 19771: 
Weige lt & Wirnitzerl Il983l) then provides accurate phases at 
high spatial frequencies. The narrow-band image reconstruc- 
tion was carried out using the code by JanBen ( 2003), which we 
employ aft er minor mo difications. This c ode implements the 
method of Keller & von der Luhel dl992l) . The instantaneous 
Optical Transfer Function is obtained from the instantaneous 
broad-band images assuming that the broad-band speckle re- 
constructed image is the true scenery. Its Fourier transform is 
denoted by Ob, where the hat indicates that we deal with an es- 
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Fig. 1. Comparison between a Stokes / profile from our flat 
field data (open circles) and a reference solar spectrum be- 
fore (dotted line) and after (solid line) convolution with the FPI 
band-pass function (FWHM 44 mA). The wavelength A is re- 
ferred to 6300 A, and the Stokes / profiles are normalized to 
the continuum intensity I c . 

timate. Under this assumption, one can write down an equation 
that yields the Fourier transform of the reconstructed narrow- 
band image 0„ as a function of known quantities, namely, Ob, 
and the Fourier transforms of the individual broad-band hj and 
narrow-band I„j images, 



2j j Injlh ; „ 



(1) 



The sums comprise all the images taken for a given wavelength, 
H represents a noise filter, and the superscript * denotes com- 
plex conjugate. All the symbols in Eq. Q are functions of the 
spatial frequency (not explicitly included). Equation Q pro- 
vides the base for the reconstruction technique. It derives, apart 
from the noise filter H, from a least-squares fit that minimizes 
the merit function 



(2) 



with respect to O n . The instantaneous Optical Transfer 
Function S j has to be replaced by 



hi. 

6 b 



(3) 



fl992 i. Krie g et alJ £999), 



or 



see Keller & von der Liihe 
Koschins kv et all d2001lh 

The noise filter H is usually an optimum filter calculated ac- 
cording to the signal (e.g., Kri eg et alJ Il999h . However, the sig- 
nal depends on the wavelength within the spectral line, which 
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Fig. 2. Filters and power spectra involved in the process of 
image restoration, a) Solid line: noise filter applied to the se- 
ries of individual frames that give rise to each restored image. 
Dashed line: optical transfer function for a long exposure im- 
age with Of.' 5 FWHM seeing, b) Dotted line: power spectrum of 
the broad-band image. Dash-dotted line: power spectrum of the 
narrow-band continuum image processed with the noise filter. 
Solid line: the narrow-band continuum image after smoothing. 
Dashed line: power spectrum of the narrow-band continuum 
image considering long exposure seeing. Frequencies are given 
in 1/arcsec. The corresponding periods are also shown on the 
upper bound of the plot. 



would produce a noise filter varying along the wavelength scan 
and, consequently, a spatial resolution different for different 
wavelengths. In this study we need to assure the same res- 
olution for all wavelengths since the reconstructed images at 
different wavelengths are combined to form individual spectra. 
For this reason we use a single noise filter which equals one un- 
til the spatial frequency corresponding to 075, and then it drops 
to zero from this point on (see Fig. The 075 cutoff was se- 
lected because it is not far from the values provided by the op- 
timum filter, and it still renders a good angular resolution. To 
further suppress the noise, the final narrow-band images were 
smoothed with a 5 x 5 pixels boxcar. Figure|5]shows azimuthal 
averages of various power spectra that take part in the pro- 
cess of speckle restoration. Note, in particular, how the broad- 
band image has power up to a spatial scale equivalent to 0725, 
whereas the corresponding power spectrum of the narrow -band 
image drops off beyond 075. 
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Fig. 3. Example of reconstruction, a) Broad-band image (spatial resolution ~ 0725). b) narrow-band continuum image (spatial 
resolution ~ (X'5). c) Intensity at the core of Fei /16302.5 A. The axes are in arcsec. 



Figure [3] shows a reconstruction. The axes are in arcsec so 
the FOV corresponds to 14" x 23". Figure [3^ contains the 
broad-band image. The contrast of this image turns out to be 
12.5% (by definition, the contrast is the standard deviation 
of the intensity in the image normalized to the mean inten- 
sity). Figure [3J> shows the continuum intensity of the narrow- 
band reconstruction. In this case the contrast drops to some 
7%. Figure |5J; is the intensity at the core of Fei /16302.5 A. 
Incidentally, the latter image provides yet another evidence 
supporting that the observations correspond to a very quiet Sun 
region. Network concentration s show up as co nspicuous bright 
points at line core images fe.g.. lSheelevll967l) . 

3.2. Stokes I and Stokes V profiles 

The restoration process is carried out independently for the two 
images corresponding to the right and left circularly polarized 
beams S i and S 2- Once the restoration is completed, we com- 
pute the Stokes / and Stokes V profiles 2 by adding and sub- 
tracting the two images at each wavelength, 



V oc S i - $2, 
/ oc S l + 52. 



(4) 



This task demands a careful superposition of the beams, which 
we carry out after sub-pixel interpolation. An example of a 
Stokes / profile is included in Fig.^ Several Stokes V profiles 
are shown in Fig.0] 



2 The term Stokes profiles denotes the variation with wavelength of 
the four Stokes parameters. In the present work we only deal with the 
intensity / and the circular polarization V. 
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Fig. 4. Set of representative Stokes V profiles of 
Fei /16302.5 A. Despite the noise, deviations from the 
anti-symmetric shape are evident. (The large plus signs point 
out the origin of abscissae and ordinates.) We do not know 
which fraction of these asymmetries is real and which results 
from contamination with Stokes /. Wavelengths are in mA off 
line core Aq. The polarization signals have been normalized to 
the continuum intensity I c . The flux density assigned to each 
profile is included for reference. 



3.3. Magnetograms and their calibration 

Magnetic flux densities are computed starting from the the 
well-known magnetograph equation, which relates the circular 
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polarization V and the longitudinal component of the magnetic 
field B (e. g.. lUnnol 19561 lLandi Pegl' Innocentil 1 992t) . 

V(A) = C(A) B, (5) 

with 

C(A) = k 8L A 2 d -^. (6) 

The symbols A, Aq and gi stand for the wavelength, the cen- 
tral wavelength of the line, and the effective Lande factor, re- 
spectively. When Stokes V and / have the same units, the con- 
stant k is equal to - 4.67 • 10~ 13 A -1 G We apply an inde- 
pendent calibration to each position of the field of view. First, 
the derivative of Stokes / required to evaluate the calibration 
constant C(A) is computed numerically. Then, we select the 
wavelengths of the two extrema of the derivative, plus the two 
wavelengths immediately adjacent to each one of them. The 
selection renders six independent wavelength positions with 
C(A) + where Eq. Q could be applied to estimate B. Each 
wavelength gives a different value, so we chose a best estimate 
solving Eq. (|5j by lest-squares. This procedure defines our es- 
timate of the flux density B e fj , 

XjViAjWj) 



AB eff ^ 17 G for Fe 716302.5 A. 



(10) 



B 



eff 



ZjCHAj) 



(7) 



The sum considers the six wavelengths Aj whose selection is 
described above (see also the square symbols in Fig. IB.2I . 
Obviously, if the observed Stokes I and V profiles follow 
Eq. l|5} then B e ff = B. In practice this is not the case, 
and an extensive literature discusses various bias resulting 
from the brea k down of the magnetograph equation (refer to, 
e.g..lJefferies & Micke\ll99lllKeller et alJl994HGraham et al) 
2002). Here we need to consider the main bias arising when the 
magnetic structure is much smaller than the resolution element, 
mostly filled by unmagnetized plasma. In this case, 



Beff * aB, 



a being the fraction of resolution element occupied by the oth- 
erwise uniform magnetic structure. According to Eq. (|8j, # e ff 
turns out to estimate the magnetic flux across the resolution el- 
ement divided by the area of the resolution element (i.e., the 
magnetic flux density). 

Two main sources of noise limit the precision of the flux 
density derived from Eq. @. The measure is affected by the 
random noise of the Stokes V spectra. We estimate its influence 
applying the law of propagation of errors ("e.g.. lMartinll97lh to 
Eq. Q, i.e., 



AB le * J] [dB en /dV(Aj)fAV(Aj) 2 



AV 2 /Yc 2 (Aj), 



(9) 



where one assumes that the noise of the different wavelengths 
AV(Aj) is independent but has the same variance AV 2 . We have 
applied the previous equation to all individual profiles. The 
mean values of the errors are 

AB eff a 23 G for Fe / 16301 .5 A, 



The AV required to evaluate Eq. (|9jl was set to 5 x 10~ 3 I c , I c be- 
ing the continuum intensity. This figure was estimated from the 
standard deviation of V as a function of A when Z? e ff becomes 
zero. Explicitly, the standard deviation of all observed spectra 
was represented versus B e ff ■ Such scatter plot gives a non-zero 
Stokes V standard deviation for B e ff = 0. This value is used for 
AV. 

A second source of error affecting B e ff is due to the con- 
tamination of the polarization signals with intensity. Following 
the procedure described in Sect. 13.21 spatially restored spectra 
proportional to I + V and / - V are obtained for each point on 
the solar surface, 



Si = /i (/+V)/2, 
S 2 = / 2 (/-V)/2. 



(11) 



The circular polarization is estimated by subtracting them out 
(Eq.ffl), 



V = Si-S 2 



(12) 



This would yield the true solar V signal if /i - fi = 1. 
However, due to the uncertainties of the reduction procedure 
(e.g., normalization to the continuum intensity, insufficient flat- 
fielding, non-linearities of the cameras, instrumental polariza- 
tion induced by the telescope, etc), f\ + fi and V is contami- 
nated by the intensity profile I, 



- h-h 
V ~ v + - — — /. 



(13) 



Since V <*c /, the crosstalk with intensity seriously threatens the 
polarimetric accuracy of the measurements. Equation Q was 
devised to automatically correct for this crosstalk, at least to 
first order. One can estimate the residual error due to crosstalk 
with intensity as follows. Equations Q and Jl 3i lead to 



(8) Afl eff = B : 



' eff ' 



B 



eft'O 



-AB, 



(14) 



with Befjo the flux density to be retrieved if there were no 
crosstalk, and AB characterizing the bias, 



AB 



ZjCHAj) 



(15) 



Note that AB only depends on the intensity. We have estimated 
its value using a mean quiet Sun profile smoothed and sam- 
pled to mimic the observational procedures: AB ~ 150 G for 
Fei 16301.5 A and ~ 120 G for Fei 16302.5 A. Keeping 
in mind that the crosstalk ( f\ - h)/2 has to be smaller than 
a few per cent (otherwise the contamination would exceed 
the real Stokes V signals), the bias induced by the crosstalk 
with intensity AB (f\ - h)/2 is at most a few G. This fig- 
ure is negligibly small compared to the error produced by ran- 
dom noise (Eq. (I10» . Should we had used a different proce- 
dure to measure B e ff, the bias would be unbearable. For exam- 
ple, using only one wing of the line AB ~ 3.5 x 10 3 G for 
Fe 1 16301 .5 A, which renders AB e ff ^ 70 G for a 2% crosstalk 
(L/i _ /i]/2 = 2 • 10~ 2 ). The huge difference with respect to 
our procedure is due to the cancellation of positive and negative 
contributions in Eq. ( fTBl when the two wings are considered. 
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Fig. 5. Magnetograms from Fei ,16302.5 A (left) and 
Fe i /16301.5 A (right). The signals below the noise level are set 
to zero. Dark and light represent different polarities, as coded 
in the vertical bar. Tick marks are separated by 1 arcsec. 

3.4. Velocities 

A brief discussion on the velocities associated with the mag- 
netic signals is included in Sect. 14. 31 These velocities are com- 
puted as the wavelength of the minimum of a parabola that fits 
the core of the Stokes / profiles. The zero of the velocity scale 
is set by the mean velocity across the full FOV, so it is affected 
by the convective blueshift. 

3.5. Time series 

We analyze the time variation of the magnetic signals (Sect.[5J. 
In order to produce a movie from the individual speckle re- 
constructed magnetograms, the individual snapshots of the se- 
ries were co-aligned allowing for a global shift among them. 
Such displacement is computed by minimizing the difference 
between the successive broad-band images of the time series. 
The mean broad-band image of the series constructed in this 
way has a contrast of the order of 7%. This figure has to be 
compared to the contrast of the best image in the series, which 
is some 12.5% (Sect.l3~Tl. 

4. Results 

4.1. Mean flux density of the magnetograms 

Figure[5]presents two simultaneous magnetograms taken in the 
two iron lines. (They correspond to the snapshot of the time se- 
ries having the best angular resolution, which is the one that we 
analyze in detail.) They exhibit a salt and pepper pattern with 
patches of opposite polarity in close contact. Approximately 
60 % of the FOV contains polarimetric signal above noise in 
one of the spectral lines. This large area coverage, together 
with the fair magnetic sensitivity of the measurements, yields 
more magnetic flux than any previous observation of the pho- 
tospheric IN magnetic fields (Sect. I6.2l i. In order to quantify 



the amount of flux, we have evaluated the mean unsigned flux 
density in the observed region \B t 



effl, 



Jeffl 



N ^ 



B' 



effl- 



(16) 



where the sum includes only those pixels with signal above 
noise in one of the spectral lines, and N represents the total 
number of pixels in the FOV. (|B e ffl is the average |B e ff| over the 
FOV once signals below noise are set to zero.) The result of the 
computation yields 



\B, 



effl 



21 G, for Fe 1/16301.5 A; 
17 G, for Fei ,16302.5 A. 



(17) 



These mean flux densities are weakly affected by the noise of 
the individual B' s , since |fi e ff I results from the average of thou- 
sands of individual pixels. The fact that we average absolute 
values and the use of a threshold do not modify this argument. 
We carried out Monte-Carlo simulations that show the robust- 
ness of the estimates d!7i . Using the full magnetograms (in- 
cluding effective flux densities below the noise), we added ran- 
dom Gaussian noise according to our estimate of the noise level 
(Eq. dlOli). The application of the definition dl6> to different re- 
alizations of magnetograms with mock noise renders the same 
mean flux densities within a fraction of one G. 

We have also computed the mean signed magnetic flux den- 
sity for the magnetograms in Fig. |5]( with a definition similar to 
the unsigned flux Eq. d 1 6t > but using B'~ instead of l*Bg ff D- It 
turns out to be, 



B, 



eff 



+2 G, 
+3 G, 



for Fe i ,16301.5 A; 
for Fe i ,16302.5 A. 



(18) 



Note that B e ff <k |B e ff I, which probably constrains the physical 
mechanisms responsible for the observed magnetic structure 
(they have to create large unsigned magnetic fluxes and, simul- 
taneously, low signed fluxes). This constraint is probably even 
more severe than the values inferred from combining Eqs. d!7i 
and (II 81 . since we cannot discard systematic errors of a few G 
affecting B e g (see, e.g., Sect. 13. 3> . 

4.2. Magnetic field strengths and filling factors 

The polarization signals obtained from the two spectral lines 
are correlated (see Fig.|5jl. However, the magnitude of the effec- 
tive flux density of Fe 1/16301.5 A, B e ff(6301), is systematically 
larger than the effective field derived using Fei /16302.5 A, 
Z? e ff(6302). The difference is inferred from mean fluxes (e.g., 
from Eq. ilH). since the noise prevents conclusions based on 
individual measurements 3 . We calculate the mean excess of 
B e if(6301) with respect to B e ff(6302) by means of a lest-squares 
fit that accounts for the errors of both B s s (6301) and B e ff(6302). 



3 The error of the ratio based on the observations in a single pixel 
is A[B ef r(6301)/fi e ff(6302)] =s 30G/|£ eff (6302)|, which follows from 
the law of propagation of errors, Eq. <10> . and the ratio of effective 
fluxes that we will infer (Eq. {T3J). Since |B cff (6302)| < 150 G, then 
A[B cff (6301)/B cff (6302)] > 0.2, which does not suffice to detect the 
ratio 1 1 91 . One needs to average over many pixels to distinguish the 
ratio from unity. 
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We only consider flux densities above noise. Our best estimate 
is 

fl e ff(6301)/fl eff (6302) 1.25 ± 0.14. (19) 

The error bar corresponds to the standard 68% confidence level, 
and it is worked out in Appendix |A] where we also reject 
the hypothesis B e ff(6301) = B e ff(6302) with high confidence. 
The 25% systematic difference can be easily interpreted if the 
intrinsic magnetic field strengths of the structures that give 
rise to the signals in the magnetograms exceeds 1 kG. Should 
the magnetic fields be intrinsically weak (say, a few hundred 
G), the magnetograph equation holds since the basic condi- 
tions for the approximation to be valid are satisfied. First, the 
Zeeman splittings of both Fei /(6301.5 A and Fei /16302.5 A 
are smaller than the Doppler width of the lines. Second, the 
magnetic fields are dynamically weak, and they cannot mod- 
ify the thermodynamic conditions with respect to those in the 
un-magnetized atmosphere. Even if one does not spatially re- 
solve the magnetic structures, the calibration constants derived 
from the observed Stokes / are valid, Eq. © applies, and 
B ef f (6301) /B eff (6302) = 1. Consequently, the fact that the 
observed ratio d!9i differs from one implies the existence of 
strong fields. These qualitative arguments, in the spirit of the 
lin e-ratio method o f Stenflo ( 1973), we re alre ady put forward 
bv lSocas-Navarro & Sanchez Almeida] d2002t) . The arguments 
can be made quantitative computing the ratio of effective field 
strengths in different model atmospheres whose intrinsic field 
strengths are known. We have carried out such calibration, and 
the details are given in AppendixlBl The ratio of effective mag- 
netic field strengths versus the intrinsic field strength is com- 
puted in a variety of atmospheres (Fig. IB. it . Two main results 
arise from this calibration: 

- The observed ratio dl9> indicates the existence of mag- 
netic fields larger than 1 kG in the photospheric lay- 
ers where the observed polarization is formed. More 
precisely 900 G < B < 1600 G when 
1.1 < B e// (6301)/B e/ /(6302) < 1.4. 

- The fact that the ratio d!9i corresponds to kG fields is al- 
most insensitive to the details of the model atmosphere. 
This property was expected according to the qualitative ar- 
guments given above. 

Note that these conclusions do not discard the existence of 
sub-kG magnetic field strengths in our data. Despite the fact 
that the characteristic field strength seems to be kG, the scatter 
among the individual ratios is very large. Consequently, there 
are many individual points whose ratio can be unity or smaller 
which, according to the calibration above, implies sub-kG field 
strengths. 

Structures with intrinsic kG magnetic field strength show- 
ing 20 G flux density have to occupy only a small fraction of 
the solar surface. The simple 2-component model allows to 
estimate the area coverage or filling factor. The filling factor 
is just the ratio between effective and intrinsic field strengths 
(Eq. ©). Iffi- lkGthen 

a ~ B e «/B ~ 20 G/1000 G ~ 0.02. (20) 
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Only 2% of the solar surface produces the observed signal. We 
observe them to cover 60% of the FOV because of the limited 
spatial resolution of the observations. One can also employ the 
same order-of-magnitude calculation to estimate the size / of 
the magnetic concentrations. If M magnetic concentrations oc- 
cupy a resolution element of size L, then 

a ~ (l/L) 2 M, (21) 

which, together with Eq. (|8j, renders 

r75km forM=l; 
/ ~ Ly/B e if/(MB) ~ | 25 km forM=10; (22) 
llOkm forM=50. 

We use L = 075 = 363 km, B - I kG, and B eff = 40 G, the 
latter being the typical value of the signals above noise in the 
magnetograms of Fig. (|5}- 

4.3. Spatial distribution and brightness of elements 
with magnetic signal 

Figure [6] shows the speckle reconstructed broad-band image 
and, overlaid to it, the Fei /16302.5 A magnetogram. Most of 
the magnetic fields are located in intergran ular lanes. Th i s ten- 
dency has been observed be fore dLin & Rimmelelfl999t rLites 
2002| ISocas-Navarrol Eo03l) and is exp ected from numerical 
simul a tions of magneto-convectio n re.g.. lWeissll978tlCattaned 
1999; V ogler & Schiisslerl 12003!) . However, some magnetic 
field signals are also found in granules , in a greement with the 
observ ations of IStoloe & Kneerl J2000l) and Koschi nskv et alJ 
(2001). This bias of the magnetic signals towards intergran- 
ules is also evident from the inspection of Fig. which con- 
tains histograms of the distributions of continuum intensities 
and velocities. The figure displays both, pixels with signals 
above noise (the solid lines), and below noise (the dotted lines). 
Some 65% of the magnetic signals of Fe i /16302.5 A are darker 
than the mean intensity, which indicates association with inter- 
granular lanes. On the other hand, 60 % of these pixels present 
redshifts, which correspond to downflows in the solar atmo- 
sphere and therefore trace intergranules. The histograms based 
on Fe i /16301.5 A present similar trends. 

The strong signals in the magnetograms show a pattern 
whose characteristic size corresponds to the mesogranular con- 
vective cells, i.e., with a spatial scale between 5 and 10 arcsec 
(Fig.lHJ. Mesogranulation has been detect ed in velocity (e.g., 
[November et al.ll98ll) and intensity (e.g.. lDeubnerll 989). and 
it is also recover ed from numerical s imulations of magneto- 
convection (e.g., ICattaneo et al.l 120011) . However, to the au- 
thors' best knowledge, this is the first clear detection in polar- 
ization. (For another observation suggestive of mesogranula- 
tion in polarized light, see Truiillo Bueno 2003, Fig. 4.) Figure 
[8] left, portraits the magnetogram of the region except that 
only fluxes above 50 G have been represented. One can guess 
a cellular pattern with a scale larger than the granules but still 
smaller than the supergranulation (the latter with a scale larger 
than the FOV). This pattern is long-lived since it shows up in 
the mean magnetogram of our 17 min time sequence (Fig. [8] 
right; see also Sect. |5} Our finding should be regarded as an 
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Fig. 6. Speckle reconstructed broad-band image overlaid 
with the magnetogram of Fei /16302.5 A with contours at 
Bett = ±30, ±50, ±70, and ±90 G. The solid and dotted con- 
tours indicate opposite polarities. The distance between tick- 
marks is 1 arcsec, as indicated by the labels of the axes. The 
arrow points out an isolated bright point without associated po- 
larization signal. 

additional observational proof for the existence of the meso- 
granular convective scale. 

We have found magnetic structures that typically harbor 
kG magnetic fields. According to the common wisdom, they 
should be bright in broad-band images. The presence of a kG 
field reduces the density of the magnetized plasma and thus 
lowers the opacity with respect to the ambient atmosphere. 
One sees deeper through the magnetized plasma which (usu- 
ally) means observing hotter layers that p roduce more light 
(often called hot-wall effect; Spruit 1976). Except for a few 
cases (e.g., bright points at [T.'5, 22"], [12", 1"], or [8", 10"] 
in Fig.[6ji, we do not see bright points associated with the mag- 
netic signals. Moreover, the intergranules having magnetic sig- 
nals do not seem to be brighter than those without them (see 
the histograms in Fig.|5Jl. This seeming inconsistency between 
the expected but unobserved brightness may not be real. There 
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Fig. 7. Histograms of continuum intensities (a) and veloci- 
ties (b) for pixels with magnetic fields above the noise level 
(the solid lines) and below the noise level (the dotted lines). 
The intensities refer to the mean intensity, and the average 
velocity over the FOV is set to zero. We adopt the conven- 
tion that redshifts are positive. The histograms correspond to 
Fe i /16302.5 A. 

are several ways to reconcile the observations with the current 
paradigm, for example, 

- the sizes of the individual kG magnetic concentrations (see 
Eq. i22\ ) are small compared with the sizes of intergranular 
lanes. Even very bright but small magnetic structures may 
lead to dark signals when observed with the kind of (good 
but still in sufficient) angular res olution of the broad-band 
image (see Title & Berger 1996). 

- If the temperature of the intergranules is low enough, the 
expected reduction of opacity associated with the presence 
of kG magnetic fields may not produce a bright struc- 
ture. The depression of the observed layers by a few tens 
of km may not be enough to reach temperatures larger 
than thos e of the mean u n-depressed photosphere (see, e.g., 
IStein & Nordlundll998l Fig. 14). 

- Due to a yet unknown physical process, there may be a ten- 
dency for the magnetic structures to stay in the darkest in- 
tergranular lanes (e.g., they accumulate in the strong down- 
flows that tend to be particularly dark.) 

Note, finally, that the signals in Fig.|6]cannot be misiden- 
tified network magnetic concentrations, since they occur at a 
characteristic scale much smaller than the supergranulation. 

4.4. Stokes V asymmetries 

The Stokes V profiles observed in quiet network and inter- 
network regions show large deviations from the anti-symmetric 
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Fig. 8. Left: magnetogram of the region showing only strong 
signals, explicitly, flux densities above 50 G. The rest is set 
to zero. Note the regular pattern with a size similar to the 5- 
10"scale of the mesogranulation (tick-marks correspond to 1 
arcsec). Right: mean magnetogram of the time series, which 
still shows a pattern similar to that of the single snapshot. In 
this case we set to zero those flux densities below 25 G . 
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Fig. 9. Distribution of broad-band intensities in intergranules 
with and without magnetic field. Intergranules are selected as 
those points with intensity smaller that 90% of the mean inten- 
sity. The two types of line distinguish the histograms obtained 
for points with magnetic field above the noise (the solid line) 
and those with no detected signals (the dashed line). There is 
no obvious difference. The histograms provide the percentage 
of points in each bin. 



shape (|S^nchez Alme da et al,||l996[ |Grossmann-Doerth e{ all 
1996[ ISigwarth et al.lll999t ISanchez Almeida & LiteslEToOO) 



These so-called asymmetries carry meaningful informa- 
tion on the structure of the magnetic atmosphere (e.g., 
ISanchez Almeida! 19981 and references therein). It would have 
been desirable studying the asymmetries of our profiles in some 
detail. However, the Stokes I to V contamination described in 



Sect. 13.31 Eq. Jl 3i . masks the real asymmetries in a significant 
way, and we refrain from analyzing them. Once this caveat has 
been explicitly pointed out, we also would like to mention that 
asymmetric Stokes V profiles are found everywhere in the FOV. 
In particular, we observe Stokes V profiles with only one lobe 
and profiles with three lobes (see Fig.@J. 

5. Temporal evolution 

The duration (17 min) and cadence (50 s) of the series of mag- 
netograms are far from optimum to allow a detailed study of 
the temporal evolution of the structures. In addition, the po- 
larization signals are close to the noise level, and the seeing 
conditions vary along the sequence. These two factors induce 
spurious time variations which complicate the analysis. Despite 
these drawbacks, several conclusions on the time evolution 
of the magnetic signals can safely be drawn. We extracted 
them from the inspection of the sequence of magnetograms re- 
centered as described in Sect. 13.51 

Many individual features maintain their identity in succes- 
sive time steps of the series (see Fig. ^TJ| where time increases 
from top to bottom and from left to right). In fact, some of the 
signals live longer than the total time span, as it is proven by 
Fig. [8] where the strongest signals in the mean magnetogram 
of the series resemble those of any individual snapshot. This 
mean magnetogram has been computed as a weighted average 
of the individual magnetograms of the series, with the weights 
given by the inverse of the squared noise of the magnetograms 
(Eg. HOt . In order to quantify the fraction of unsigned flux that 
survives more than the time-span, we have computed the un- 
signed flux density for signals above noise level in the mean 
magnetograms. They turn out to be, 



\B. 



effl 



13 G, for Fe 1/16301.5 A; 
9G, for Fe i ^6302.5 A, 



(23) 



They represent 60% of the flux density in the best individual 
magnetogram (Eq. ilH ). It is tempting to interpret this figure 
as the fraction of IN magnetic structures that live longer than 17 
min. However, the interpretation is not so straightforward since 
we detect polarization signals rather than magnetic structures. 
Seemingly long-lived signals may not necessarily correspond 
to stable magnetic structures. Numerical simulations of quiet 
Sun magneto-convection show persistent downdrafts that con- 
tinuously gather magnetized plasma in spe cific locations (e.g. 
ICattaned ll999: Emonet & Cattaneo l200ll) . The plasma sinks 
down and disappears along the downdrafts, and it is continu- 
ously replaced by new plasma. The system of persistent down- 
drafts outlines mesogranular cells, which live much longer than 
the individual magnetic structures. This theoretical scenario 
is consistent with the observed magnetograms and, in partic- 
ular, with the existence of a long lasting polarization pattern 
that we have already associated with the mesogranulation (see 
Sect. 14.31 1. Consequently, the existence of signals in the mean 
magnetogram is interpreted here as a proof of the temporal co- 
herence of the mesogranular pattern, rather than the persistence 
of individual magnetic concentrations. 

On top of the rather stable pattern described above, large 
variations of the magnetic signals occur between snapshots 
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Fig. 10. Temporal evolution. Magnetograms overlapped to intensity images as in Fig. ©but for a sub-field of 7" x7". Time 
increases from top to bottom and from left to right. The different images are consecutive snapshots of the series separated by 50 s 
(see the labels on the images). The arrow points out a specific region dominated by a magnetic concentration of negative polarity 
that interacts with positive polarity patches and, possibly, eats them away. 
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(Fig. II Oi . Some polarization signals migrate following the 
granular motions to remain in the intergranular spaces. There 
are places where opposite polarities approach each other and 
seem to (partly) cancel out; the arrow in Fig. ^3 indicates a 
negative polarity patch (dotted line) that interacts, and possi- 
ble eats away, neighboring positive polarity contours (compare 
the snapshots at m s and 4™ 10 s ). More often signals grow 
and fade with no obvious interference with patches of oppo- 
site polarity. As we pointed out above, these variations only 
trace changes in the polarization signals, which do not neces- 
sarily imply the evolution of the magnetized plasma. Actually, 
one may have large variations of polarization without substan- 
tial changes of the underlying magnetic structure. For example, 
ISanchez Almeidal |2000) describes how the polarization of kG 
magnetic concentrations can vanish upon a small increase of 
the field strength. This effect is mentioned here because it may 
be responsible for the lack of magnetic signal associated to the 
bright point indicated by an arrow in Fig. [6] Subsequent mag- 
netograms of the series show that this bright point is co-spatial 
with polarization signals. 

The mean unsigned flux density is only slightly different 
for the different snapshots of the series. The standard deviation 
of these changes is some 10% for Fei /16302.5 A, whereas it 
reaches 20% for Fei /16301.5 A. There is a tendency for the 
snapshots of largest continuum contrast to have the largest flux 
densities. 

6. Consistency of the results 

The unsigned flux densities that we find are larger than the val- 
ues found hitherto. The purpose of this section is to spell out 
arguments that support the reliability of our results. The mag- 
netograms have self-consistency since independent parts of the 
data set provide similar results. They are also consistent with 
previous observations having different sensitivities and angu- 
lar resolutions, once these differences are properly taken into 
account. Finally, the signals meet several theoretical prejudices 
which are independent of the observation. In particular, the sig- 
nals appear in the intergranular lanes and show a conspicuous 
and stable mesogranular pattern. 

6.1. Self-consistency 

We have a time sequence of magnetograms taken in two differ- 
ent spectral lines. Magnetograms differing in time or spectral 
line are independent, and they can be analyzed independently 
to check the self-consistency of various results. We have iden- 
tified the following properties of the data set that support the 
internal agreement between independent subsets, 

- the simultaneous magnetograms obtained from the two 
lines are very similar (Fig.FJJ. 

- Many patches in the magnetograms are larger than the res- 
olution of the observation, implying a spatial coherence of 
the signals difficult to ascribe to noise. 

- The signals are co-spatial with the dark intergranular lanes 
(Sect. l4~3l . Again, this relationship is difficult to explain as 
due to noise. One might suspect that the misalignment of 



the two circularly polarized beams subtracted to determine 
Stokes V creates false polarization signals. This contami- 
nation would be maximum where the spat ial gradie nts in 
intensity or velocity are maximum (e.g.. lLitesl ll987). The 
maximum gradients occur in the transitions between gran- 
ules and intergranules, which is not where the detected sig- 
nals are. We have checked that the scatter plots of polariza- 
tion signals versus spatial gradients show no obvious cor- 
relation. 

- Many polarization patches maintain their identity in sev- 
eral successive magnetograms of the series. In particular, 
the mean magnetogram retains some 60 % of the peak flux 
density (Eq. {23}), a fact once again difficult to explain as 
due to noise. Should the signals be noise, the mean magne- 
togram would show signals much smaller than that of indi- 
vidual frames (i.e., only ~ 1 / V20 ~ 20 % of the individual 
signals). 

6.2. Comparison with other measurements 

Our observations differ in angular resolution and sensitivity 
from previous measurements. These factors have to be taken 
into account before they can be properly compared. We have to 
cut down the resolution of the magnetogram and then change 
the sensitivity of the measurement. 

The speckle reconstructed images can be regarded as the 
true solar image down to the cutoff frequency. In the case of 
our magnetograms, the cutoff corresponds to a period of 075 
(see Sect. 13.11 Fig. |3. We model the effect of seeing on the 
measured images by convolving the speckle reconstructed im- 
ages with a Gaussian point spread function. This smoothing is 
carried out independently for each single wavelength. Then the 
different wavelengths are combined and processed to get mag- 
netic flux densities in the same way as the original images. The 
FWHM of the Gaussian kernel is used to quantify the seeing. 
Figure^2 sriows how the magnetograms (left) and broad-band 
intensities (right) appear under various seeing conditions 4 : 073, 
075, 1" and 2". The field of view of the section of the mag- 
netogram that we present is only 7" x 7". Even when the see- 
ing is fair (1"), only traces of the signals in the speckle re- 
constructed images are left. The unsigned flux density that sur- 
vives a given seeing is shown in Fig. El It contains the signals 
above three different levels of noise (or different sensitivities): 
100 G, 20 G and the actual noise in the magnetograms. The de- 
crease of angular resolution reduces the signals due to cancel- 
lation between close opposite polarities. For an angular resolu- 
tion corresponding to 1"(= 0.75 Mm), the mean unsigned flux 
density become of the order of 7 G. This value agrees with pre- 
vious determinations bas ed on scanning spectro-po larimeters 
that reach 1" resolution (Sanchez Almeida & Lites 2000 and 
lLitesll2002l obtain some IPG). ISanchez Almeida et alJ fc003) 
collect different values from the literature, showing a magni- 

4 One might regard as inconsistent modeling a 073 seeing when the 
cutoff frequency of the observation is only 075. Note, however, that 
the speckle reconstructed images contain all the power down to the 
cutoff, and this power is partly reduced under any seeing condition, 
including seeing with FWHM smaller than the cutoff. 
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Fig. 11. Magnetograms (left) and continuum intensities (right) 
for different seeing conditions as indicated by the insets. All 
magnetograms have the same scale, which saturates at -60 G 
and 60 G. The field of view is only 7" x 7". 

tude and trend that resem bles the solid line in Fig.^] In par- 
ticular, IWang et al.ldl995l) point out a flux density of 1 .65 G for 
a 2" resolution, which is not far from the value that remains 
in our magnetograms for this resolution (some 2.5 G; Fie. 1121 . 
On the other hand, magnetograms with an angular resolution 
similar to those presented here exist in the the literature (Keller 
ll995UKoschinskv et al.ll200lUBerger & Tiflell200ll) . However, 
they have a noise level between 50 G and 150 G, which exceeds 
the limit required to detect most of the signals in our magne- 
tograms. This lower sensitivity can easily explain why the IN 
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Fig. 12. Mean unsigned Fei /16302.5 A flux density that re- 
mains in the magnetograms when observed with the seeing 
given in the abscissa. If all signals above noise are considered, 
one is left with the solid line. If only signals above 20 G are 
considered, then the reduction is much more severe as repre- 
sented by the dashed line. Finally, when very large flux densi- 
ties are considered (signals larger than 100 G) almost no signal 
remains (the dotted line). 



signals that we detect have been missed so far. FigurefT2lshows 
that the mean flux density drops below 1 -2 G when only sig- 
nals above 100G are considered. In short, our magnetograms 
are consistent with previous observations when the limitations 
of angular resolution and sensitivity are taken into account. 

We find imbalance between the two polarities so that there 
is a net magnetic flux across the region (Eq. dl8ll). We do not 
claim that the effect is real since it is at the level of systematic 
effects that may be important (see Sect. I3.3I >. However, imbal- 
ances of simil ar magnit ude have been found in the IN by other 
authors (see Lites 2002), and consequently, the imbalance may 
be true. If so, it would be one more aspect of our observation 
being consistent with previous measurements of IN field prop- 
erties. 
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7. Conclusions 

We have obtained a time sequence of speckle reconstructed 
quiet Sun magnetograms showing magnetic signals that cover 
60% of an Intra-Network region (Fig.|6j. The signals appear as 
patches of opposite polarity often located close to each other. 
The structures are not evenly spread out, but they tend to ac- 
cumulate in the intergranular lanes. Those signals of large flux 
density trace a network whose scale is larger than the gran- 
ules and smaller that the supergranules. We associate this net- 
work with the mesogranulation, difficult to observe in velocity 
and intensity but very conspicuous in polarization (Fig.|8j. The 
mesogranular pattern lives longer than the time span of the se- 
ries (17 min). 

The effective flux density deduced from Fei ,16301.5 A is 
systematically larger than the flux density of Fei A6302.5 A 
(Eq. We interpret this difference as due to the presence of 
kG field strengths in the underlying magnetic structures. These 
field strengths saturate, differently, the polarization signals of 
the two spectral lines. The inference of kG, however, does not 
rule out that weaker field strengths are responsible for some of 
the observed signals (see Sect. I4.2i . Since the flux density is 
proportional to both the field strength and the filling factor, we 
are forced to conclude that the true magnetic concentrations 
producing the polarization cover only a small fraction of the 
surface; some 2%. 

The magnetic signals that we observe meet several theoret- 
ical prejudices: they appear in intergranular lanes, trace meso- 
granulation, vary in short timescales, and have complex topol- 
ogy with opposite polarities in contact. These properties are 
easy to interpret as the result of the interaction between gran- 
ular convection and magnetic fields, no matter whethe r the 
fields are produced by a local dvnamo JPetroyav & Sza kalv 
1993; Cattaneo 1999; Emonet & Cattaneo 2001), a global dy- 
namo (iStein & Nor dlund 2002), or they repre sent re-processed 
materials left by old a ctive regions (e.g., ISchusslerl 120031 
IVogler & Schusslerl2003l) . There are other observational prop- 
erties which do not accommodate so easily within the existing 
paradigms. The field strengths are typically kG, so a process 
has to concentrate the fields to this level. The convective col- 
lapse ( Park erll978tlSpruitll979i) . devised for network magnetic 
concentrations, may have difficulty to concentra te the low mag - 
netic flux features that we detect (see Sanchez Almeida 2001). 
Yet another observational constraint difficult to satisfy is the 
large difference between the signed and unsigned flux densities. 
Whatever physical mechanism generates the observed fields, it 
has to produce a distribution of magnetic field whose mean flux 
is 10 times smaller than the local fluctuations of flux. 

ISanchez Almeida et alJ J2003I) compile values for the mean 
unsigned flux density of the IN fields measured by different 
authors. They show a tendency to increase as the angular res- 
olution improves, but none of them exceeds 10 G. Our magne- 
tograms contain more magnetic flux than the values reported 
hitherto. We detect more polarization signals probably due to 
the unique combination of high spatial resolution and good sen- 
sitivity. The magnetic flux detected in the IN critically depends 
on the polarimetric sensitivity and the angular resolution. If 
these factors are properly taken into account, our estimates of 



unsigned flux density are compatible with the previous mea- 
surement. 

We find evidence that our angular resolution and sensi- 
tivity do not suffice to detect all the magnetic flux exist- 
ing in the quiet Sun. The flux density in the magnetograms 
drops off upon artificial reduction of polarimetric sensitiv- 
ity or spatial resolution (Fig. \121 . The extrapolation of this 
trend predicts a notable increase of signals upon improve- 
ment of the observational limitations. On the other hand, 
we have concluded that the observed signals are tracing 
kG magnetic fields. There are sensible theoretical arguments 
pointing out that the spectra l lines used to measure are bi- 
ased towards these kG fields dSanchez Almeida & Litesl f2000 ; 
ISocas-Navarro & Sanchez Almeidal 120031) . In other words, 
magnetic concentrations of sub-kG field strength could have 
been easily overlooked by our magnetograms. These and other 
reasons suggest that the IN flux detected so far should be re- 
garded as a lower limit to the true flux. Thus any effort to im- 
prove the sensitivity and/or angular resolution is likely to be 
rewarded with new magnetic flux. 
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Appendix A: Ratio B e ff(6301)/B e ff(6302) 

We want to estimate the typical ratio of effective magnetic 
fluxes taking into account that (a) the measured magnetic fluxes 
of the two lines have errors, and (b) these errors are not negli- 
gible. We invoke several results f rom statistics, which can be 
found in standard textbooks (e.g., Martin Il97lh . For the sake 
of conciseness, and only in this Appendix, the following nota- 
tion is employed 

Bu = B e ff(6301) in the i - th pixel, 

B 2i = B e ff(6302) in the i - th pixel. (A. 1) 

The fluxes observed in each pixel have true values, B° u and 
Z?",, plus noise, ABu and AB 2 ,, 

B u = B° u + AB U , 

B 2i = B% + AB 2i . (A.2) 

We assume the noise of the two fluxes to be independent, since 
they come from separate data. We also assume a linear relation- 
ship between the true fluxes, 

B°u=mB° 2i , (A3) 
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being m the ratio that we want to estimate. The above expres- 
sion, combined with Eq. iA.2\ . renders a linear relationship be- 
tween the two observables By and Bji that explicitly includes 
noise, 



By - mB 2 i - mAB 2 i + ABy. 



(A.4) 



As usual, m is determined by means of a least-squares fit. We 
chose to minimize the merit function 



X 2 (m 



(By- mB 2i f 



(A.5) 



where <x 2 . and crL are the variances of the probability density 
functions describing the noise ABy and AB 2 i, 

crl = V{AB U } = V{By}, 

cr 2 2i = V{AB 2i ] = V{B 2i }, (A.6) 

whose expectation values are assumed to be zero, 

E{ABy] = E{AB 2i ) = 0. (A.7) 

(The symbol V{ } stands for the variance of a random variable, 
and it should not be confused with the Stokes V parameter.) N 
in Eq. ( IA.5> represents the number of pixels used in the esti- 
mate. This particular merit function has been selected because 
it presents several practical advantages: 

- E{By - mB 2 j} - 0, V\By - mB 2 j} = cr 2 . + ?« 2 cr 2 ., and so, in- 
voking the central limit theorem, each term of the definition 
\A.5\ is the square of a random variable distributed accord- 
ing to a n(0, 1) (i.e., a normal distribution with mean of zero 
and variance of one). Should the noise of different pixels be 
independent, the random variable X 2 follows ax 2 distribu- 
tion, which makes it easy setting confidence intervals and 
performing statistical tests (see below). 

- The ratio of flux densities derived from the fit does not de- 
pend on whether we estimate m or m _1 , i.e., on whether 
we start up from the relationship iA.'H or its inverse Z?". = 
m'B { ).. In this second case one would obtain m' - m~ l . 

- X minimizes the relative errors of the fit including both the 
errors of abscissae and ordinates. 

We carry out the non-linear X 2 minimization by brute force, 
representing X 2 versus m and selecting the extreme. Figure lATl 
shows X 2 (m) when cry and cr 2 i are the errors of the fluxes esti- 
mated according to Eq. JlOi . If these variances were exact, then 
the expected value of each term of Eq. \A.5\ is one and 



£{X 2 } = N. 



(A.8) 



When the original error estimates are used, the value of X 2 
at the minimum becomes slightly below this expected value 
jA.8> . which we interpret as an overestimate of the variances 
when using Eq. fllOi . We cure the excess decreasing all the er- 
rors by 10% to force that X 2 equals N at the minimum. Figure 
JA. II shows that such minimum of X 2 occurs when 



(A.9) 



m= 1.25 + 0.14, 



where the error bar represents the standard 68% confidence in- 
terval. In order to evaluate this confidence interval, one must 
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Fig.A.l. Merit function minimized to estimate the ratio 
Z? e ff(6301)/B e ff(6302). This ratio is represented in abscissae us- 
ing the symbol m. The merit function X 2 has been normalized 
to the number of points used in the estimate, and it has a min- 
imum around m ^ 1.25. The horizontal and vertical lines give 
the confidence intervals for two significance levels, as indicated 
by the inset. 

know the probability density function that describes the merit 
function X 2 considered as a random variable. If all the points 
of the magnetograms were independent, then X 2 would be dis- 
tributed according to a x 1 distribution with N degrees of free- 
dom, implying 



V{X 2 



2N. 



(A. 10) 



However, we smooth the original magnetograms with a 5 X 5 
pixels boxcar (Sect. 13. 1> . so the different pixels are not inde- 
pendent. This smearing reduces the degrees of freedom of X 2 
and increases its variance with respect to that in Eq. jA.101 . 
If the magnetograms are divided in groups of 5 x 5 neighbor- 
ing pixels, the central points of all the groups are independent. 
Using only these central pixels to construct the merit function, 
then the new merit function X 2 . follows a y 1 distribution with 
A725 degrees of freedom, since now the sum (IA.5> contains 
A725 summands. If all the members of each group connected 
by the smoothing were identical, then X 2 = 25 X?, and so X 2 
would have a variance 25 times larger than that in Eq. 



V{X 2 } = V{25X 2 } 



25 2 • 2N/25 = 25 • 2N. 




Equations (IA. 101 and l lA.ll> represent the two extreme limits 
corresponding to the cases where all pixels which are not in- 
dependent due to the smoothing are identical (Eq. iA.l 1» . and 
when all the pixels are independent (Eq. JA. 101 . In general 



V{X 2 } =a 2 2N, 



(A. 12) 



with a — 1 if all the pixels have independent noise and a — 5 
if pixels linked by smoothing have identical noise. In our prob- 
lem a has to lie in between these two extrema. We estimate its 
value carrying out a Monte Carlo simulation. We produce 2D 
maps of n(0,5) random variables with the size of our magne- 
tograms. These maps are smoothed with a 5 x 5 pixels boxcar, 
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so that each individual point of the smoothed maps is a n(0,l) 
random variable (i.e., the same kind of dependent random vari- 
ables that contribute to X 2 ). Then we randomly select N pixels 
of the map to compute X 2 according to the definition dA.5l >. The 
standard deviation of X 2 among several hundreds independent 
realizations of the numerical experiment renders 



spheres is 



V m (A)= \[l q (A + A B )-I q (A-A B )], 
I,„(A)= \[l q (A + A B )+I q (A-A B )\ 



(B.l) 



a ^ 2. 



(A.13) 



where I q represents the Stokes / line profile produced by the 
unmagnetized quiet Sun. The symbol A B stands for the Zeeman 
splitting of the line, 



In this intermediate case X 2 does not necessarily follow a x 2 ^ B ~ ^ L ^fi' 
distribution, however, invoking the central limit theorem, one 
can argue that the random variable defined as 



(B.2) 



Y = (X 2 -N)l(a^2N), 



(A. 14) 



follows a n(0,l). This argument, together with the value of a 
let us set confidence intervals. For example, 



given in ( IA. 13i 
the standard 68% confidence interval, corresponding to |F| < 1, 
yields the values of m included in Eq. ( IA.9b . They are computed 
as those m with X 2 (ra) - N + a s/2N (see Fig. lA.U . 

Equation JA.14I can also be used to reject the hypothesis 
m — 1 with high confidence (i.e., to reject that the effective flux 
densities of the two lines are identical). Assume that m - 1 and 
that the error estimates are correct. The high value of X 2 (l)/N 
in Fig. lA.ll mav be produced by an unusual fluctuation of the 
random variable in our particular realization. Using Eq. JA.14> . 
the probability of having a X 2 equal or larger than the value that 
we measure is smaller than 1 (T 4 , since [X 2 ( l)-N]/a [2N] 1 ^ ~ 
4. In other words, the two lines show different magnetic flux 
densities with 99.99% confidence. 

The estimate jA.9> uses those pixels in the magnetograms 
with signal above noise. As we have mentioned, era and era 
come from the values obtained from the individual estimates 
of noise described in Sect. 13.31 However, we tried other possi- 
bilities: (a) using all pixels including those with signals below 
noise, (b) using era and era constant and given by the mean val- 
ues obtained in Sect. 13. 31 (c) forcing era = era, and (d) adopting 
a = 3.5. In all these other cases m ~ 1.2 and the chance m = 1 
turns out to be improbable. 

Appendix B: Calibration of the magnetic field 
strength determination 

We want to calibrate the relationship between the ratio of flux 
densities derived from Fei /16301.5 A and Fei /16302.5 A, 
fi e //(6301)/fi e/ /(6302), and the intrinsic field strength. First, 
Stokes / and V profiles are synthesized in a set of model atmo- 
spheres of known field strength. Second, the synthetic profiles 
are smeared out with a band-pass filter similar to the PFI used 
in the observations. Third, the profiles are sampled in wave- 
length according to the observed profiles. Finally, the definition 
of B e ff is applied (Eq. @). 

Synthetic ratios are presented in Fig. IB. II Three different 
types of model atmospheres are considered. The first one as- 
sumes a magnetic atmosphere with thermodynamic properties 
identical to that of the quiet Sun. Then, under the assumption 
of longitudinal magnetic field 5 , the polarization of such atmo- 



where the symbols are those defined for Eq. 0. Since the 
magnetized plasma may not completely fill each resolution ele- 
ment, the Stokes / and V profiles must include a fraction (I -a) 
produced by the unmagnetized plasma in the resolution ele- 
ment, 



V(i) = aV n (A), 
1(A) = al m + (1 - a)I q (A). 



(B.3) 



The solid line in Fig. lB.ll has been computed using these equa- 
tions, together with the quiet Sun prof iles I„ observ ed at disk 
center (Brault & Neckel 1987, quoted in lNeckell 19991) . We vary 
the longitudinal magnetic field strength B along the sequence, 
assuming aB to be constant (30 G; a value similar to the ef- 
fective magnetic field in our magnetograms). As expected, the 
ratio is one for weak fields and increases to 1.4 for B of the 
order of 1.5 kG. Figure lB~2l illustrates the reason of such be- 
havior. It shows the two terms of the magnetograph Eq. (|5} for 
weak fields (B = 150 G) and strong fields (B = 1500 G). When 
the field is weak then the magnetograph equation is a good ap- 
proximation and the two terms are identical. The two spectral 
lines give the true flux density and therefore the ratio becomes 
one. On the other and, the magnetograph equation is no longer 
a good representation of the Stokes V when the field is strong. 
The Stokes V profiles are much too broad as compared with 
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This assumption could be relaxed since it is not critical. 



Fig. B.l. Ratio of effective magnetic fields versus true mag- 
netic field strength in different model atmospheres. The ob- 
served ratio (Eg. 1191 ) agrees with the modeled values only if 
the intrinsic fields are larger than, say, 1 kG. This result is inde- 
pendent of the details of the model atmosphere. Different types 
of line represent different model atmospheres, as coded in the 
inset. 
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the derivative of the Stokes / profiles. (Note that Stokes / is 
almost exclusively produced by unmagnetized plasma, since 
a <K 1 in Eq. jB.3> .) This so-called saturation biases the effec- 
tive magnetic field below the real flux density. Since the effects 
increases with increasing g L , B e ft(630 l)/B e //(6302) > 1. 

In order to show that the above calibration is independent 
of the thermodynamic parameters of the atmosphere, we repeat 
the calculations using synthetic profiles from a hot model at- 
mosphere (Solanki 1986 network, with constant magnetic field, 
constant microturbulence of 1 kms -1 , and a macroturbulence 
of 1 kms 1 ). The dependence of the ratio on the field strength 
is also represented in Fig. lB.2l (the long dashed line). The dif- 
ferences with respect to the quiet Sun atmosphere are small and 
not discussed further. 

We use a third type of model atmosphere to calibrate the ra- 
tio, n amely, the model MISMAs bv lSanchez Almeida & Litesl 
(2000). These semi-empirical atmospheres are able to account 
for all Stokes profiles of Fe i ,16301 .5 A and Fe i ,16302.5 A ob- 
served in the quiet Sun at the disk center, including their asym- 
metries. The magnetic field strength of the model MISMAs 
varies with height in the atmosphere. We modify the field 
strength of the original models by changing the magnetic field 
strength at the base of the atmosphere, and then recomput- 
ing the full vertical stratification. The mean magnetic field of 
the atmosphere represented in Fig. IB. II is chosen as the mag- 
netic field strength at the height where the line core opacity of 
Fe i /16302.5 A equals one. The variation of the synthetic ratio 
versus mean magnetic field strength is qualitatively different 
than the variation found in the two previous cases, however, 
the difference is significant only when the magnetic field is 
weak. The observed ratio dl9l also requires kG magnetic field 
strengths. The reason why the ratio B e //(6301)/B e //(6302) 
does not tend to one for weak field has to do with the stray-light 
contamination of the model MISMAs (a contribution equiv- 
alent to the term [1 - a\I q in Eq. JB.3» . It has a velocity 
distribution very different from the velocities of the magne- 
tized plasma. Then the magnetized and un-magnetized atmo- 
spheres are no longer similar, which causes the break down of 
the magnetograph Eq. (|3J- Figure |B~T1 includes several kinds 
of MISMAs, from simple common ones (classes and 2) to 
atmospheres having opposite polarities in the resolution ele- 
ment (classes 6 and 7). All of them assign similar magnetic 
field strengths to the observed ratio fH?l >. 
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